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1 Introduction 



Maxwell's equations receive corrections from virtual excitations of the charged quantum fields 
(notably electrons and positrons). This leads to interesting effects [fflj: light-by- light scattering, 
photon splitting, modification of the speed of light in the presence of strong electromagnetic 
fields, and - last, but not least - pair production. 

When the heavy degrees of freedom are integrated out (in this case, the "heavy particles" 
are the electrons and positrons), an effective theory results. The corrections can be described 
by an effective interaction, the so-called quantum electrodynamic (QED) effective Lagrangian. 
The dominant effect for electromagnetic fields that vary slowly with respect to the Compton 
wavelength (frequencies u> <C 2mc 2 /h) is described by the one-loop quantum electrodynamic 
effective (so-called "Heisenberg-Euler") Lagrangian which is known to all orders in the electro- 
magnetic field |l|-|6|]. 

The Heisenberg-Euler Lagrangian AC, which constitutes a quantum correction to the 
Maxwell Lagrangian, is usually expressed as a one-dimensional proper-time integral [see e.g. 
Eq. (3.43) in H], the notation is clarified in Sec. y below]: 
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Because there are no singularities in the first quadrant of the complex plane and because Jordan's 
Lemma may be applied to the integral, it is possible to exchange the lower and upper limit of 
integration by ry and rj + ioo, respectively. 

Although the proper-time integral (|l|) can be evaluated by numerical quadrature, it is evident 
that a representation by a convergent series expansion could have certain computational as well 
as conceptual advantages. An expansion of ([I]) in terms of special functions has been given by 
us in [0,[|. Here, we present a unified series expansion which encompasses both the real and the 
imaginary part of the Lagrangian (see Sec. ||). Also, we clarify certain technical details concerning 
the derivation of our previous results |7|,[|. In Sec. |3], we discuss the "electric-magnetic duality" 
which has recently drawn much attention |],[l(|. In particular, we elucidate different kinds of 
dual invariance and pursue the question as to whether these invariances are realised by QED 
effective Lagrangians. Before we expand on these aspects, we would like to provide some general 
discussion on the general relevance of studies related to the QED effective action (|]) . 

The real part of the Lagrangian (Q) can be used to delineate dispersive phenomena such as 
photon propagation in a magnetic field (see |Qj and references therein) , photon splitting ]IT|-[l4 ] , 
vacuum birefringence and second harmonic generation [16|, and light scattering in a vacuum |j]. 
These applications have a strong relevance to particle astrophysics. The imaginary part of the 
Lagrangian (|l]) has been applied to absorption processes such as electron-positron pair creation 
and dichroism (see e.g. J4|,|7, 17-20|). The occurrence of strong fields in storage rings, pulsars, 
magnetars and high-intensity lasers also motivates our study. The advent of state-of-the-art laser 
beams and photon detectors may provide a signature of "QED's nonlinear light". The optical 
second harmonic generation (SHG) in vacuo is an interesting second-order magneto-optical effect 
that occurs if the spatial symmetry of the nonlinearity induced by the effective action (|]) is bro- 
ken e.g. by a strong static magnetic field (in a more general context, higher-harmonic generation 
by vacuum effects was discussed in fll5|,[lql ). Note that the leading contribution to SHG in vacuo 
involves a fourth-order effect in contrast to the much weaker 6th-order effect (hexagon graph) 
which gives the leading contribution to photon splitting. SHG and photon splitting in vacuo 
may be occurring close to astronomical objects such as white dwarfs, neutron stars and "mag- 
netars" (see |l4|,pn| and references therein) which have strong magnetic fields up to B ~ 10 G. 
In all cases, a detailed, realistic description of the experimental conditions and/or the involved 



astrophysical objects, especially at extreme field strengths, requires techniques for the reliable 
numerical evaluation of the QED effective action. 

A priori, the construction of a series expansion for (|l]) constitutes a complete solution of 
the problem from a theoretical point of view. However, such an expansion does not necessarily 
provide all answers: Many series expansions are known which either converge extremely slowly 
or which do not converge at all. Moreover, in the case of the Heisenberg-Euler Lagrangian, there 
is the additional problem that the terms of the series are represented by special functions which 
are in most cases defined and computed via series expansions. Again, convergence problems are 
more likely the rule than the exception. Here, we are concerned with the solution of the principal 
numerical difficulty associated with the the slow overall convergence of the series expansion 
derived in J5J|, whose terms are nonalternating in sign. The modern theory of nonlinear sequence 
transformations which begins with Wynn's epsilon algorithm |36|, |37f was developed when the 
first computers became generally available. Pade approximants, however, would not be powerful 
enough to sum our series for the Heisenberg-Euler Lagrangian derived in J7|]. 

From a mathematical point of view, there is a distant analogy between the expansion of 
(|l]) in terms of special functions and the (exact) expansion of certain quantum electrodynamic 
bound-state effects into partial waves |38|-|4l|1 . However, the mathematical entities involved in 
the present decomposition possess a far less involved mathematical structure, and it is difficult to 
associate a meaningful physical interpretation to each term in the special function representation 
of (JH). As is evident from the discussion in Sec. ||, the terms of the convergent series representation 
may be interpreted as being generated by a "partial-fraction decomposition" in distant analogy 
to the "partial-wave decomposition" in bound-state calculations. We point out in Sec. || that the 
convergence of the special function representation can be accelerated by the same technique - the 
so-called "Combined Nonlinear-Condensation Transformation" (CNCT) [Q - which was used 
successfully for the acceleration of the convergence of partial-wave decompositions in quantum 
electrodynamic bound-state calculations p0|, 41 1 . 

At the same time, we would like to mention that the integral (JlJ) , when expanded in powers 
of the electric and magnetic field strengths, represents a divergent series. The resulting divergent 
series can be used as a "model laboratory" for resummation methods |2j-25|, and related investi- 



gations [24|-27] have led to the development of asymptotically improved resummation algorithms 



which have recently received interesting applications |28|-pC|l in the highly accurate determination 
of the strong coupling constant at the Z pole and other investigations on nonperturbative effects 
in gauge theories |3l]-p4"[ . These investigations are related to the fundamental question of how 
to "make sense" of the fact that many perturbation series encountered in physics are divergent, 
and are not meant to provide efficient means for numerical evaluation of the integral ([!]) . 

2 Representation of the QED Effective Action by 
Special Functions 

The QED effective Lagrangian can be expressed as a function of the Lorentz invariants J- and 
Q which are given by 

F = 74 F, v FV = Va (B 2 - E 2 ) = % (a 2 - b 2 ) , (2) 

Q = l UF iiV (*FY v = -E B = ±ab, (3) 

where E and B are the electric and magnetic field strengths, F pi> is the field strength tensor, and 
{*FY V denotes the dual field strength tensor (*F)^ = (1/2) e^ p(J F pa . By a and b we denote 



the secular invariants, 



a 



VJ r2 + G 2 + J 7 -, 



b = \j\JT 2 + G 2 -T. (4) 

These Lorentz invariants are referred to as secular invariants because they emerge naturally as 
eigenvalues of the field strength tensor; these eigenvalues are conserved under proper Lorentz 
transformations of the field strength tensor. There are connections between the different rep- 
resentations ||]: If the relativistic invariant Q is positive, then it is possible to transform to a 
Lorentz frame in which E and B are antiparallel. In the case Q < 0, it is possible to choose 
a Lorentz frame in which E and B are parallel. Irrespective of the sign of Q we have in the 
specified frame 

a = \B\ and b = \E\ if and only if B is (anti-)parallel to E . (5) 

In any case, because a and b are positive definite, we have 

a b = \E ■ B\ > for any Lorentz frame and Q ^ , (6) 

which clarifies the sign ambiguity in (^). We give in (|2|) and (||) seemingly redundant definitions, 
but it will soon become apparent that each of the alternative "points of view" has its applications. 
The Maxwell Lagrangian is given by 

£ c , = - T = -7 4 iV FilV = l h {E 2 - B 2 ) = Va {b 2 - a 2 ) . (7) 

As it is obvious from Eq. (|l]), the correction AC to the Maxwell Lagrangian is conveniently 
written in terms of the secular invariants a and b. 

The effective action has a dispersive (real) part and an imaginary part which is associated 
with pair production, 

A£ = ReA£ + iImA/;. (8) 

In Eqs. (2)-(6) of we showed that the real part of (||) can be expressed as 

e 2 ^ 
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n=l 

coth(n7r&/a) f fnirm 2 \ ( ni:m 2 \ I ' m\m 2 \ I ' nism^ 
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n y \ ea J \ ea J \ ea J \ ea 

coth(n7ra/6) f fmrm 2 \ „ ( mrm 2 \ ( mrm 2 \ „ { rnrm 2 \ > \ . . 

d - = — k-^{ exp (^rj E '{-—) +e * p {-—) E '{—)}- <9> 

We also derived the following representation for the imaginary part [see Eq. (18) of Q]: 

^ = i?i;^(T ) exp ("^r)^ <10> 



n=l 



These results have recently been confirmed in Pj. For the special functions, we use the notation 
of Abramowitz and Stegun |43|| . Here, one might wonder why the cosine and sine integrals 
appear in an "asymmetric" form (Ci and si instead of ci and si) in the definitions of a n and d n 
in (||). The reason is that the commonly accepted definitions for the cosine and sine integrals 



are "asymmetric" in the following sense [see Eqs. (5.2.1), (5.2.2), (5.2.5), (5.2.26), and (5.2.27) 
of H]: 
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(11) 



From these formulas, it is evident that "symmetric" integrals with lower limit z and upper limit 
oo require the "asymmetric" occurrence of Ci and si. 

Because the imaginary part Im AC is generated exclusively by analytic continuation of one 
of the exponential integrals - specifically, the term 

Ei 



nirm 



eb 

in the definition of d n -, it is obvious how to write down a unified representation for both the 
real and the imaginary part. We therefore present here the unified representation for both the 
real and the imaginary part we obtained in Eq. (8) of version 2 of our preprint [44|: 
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It becomes obvious from this representation that the effective action has branch cuts along the 
positive and negative b axis as well as the positive and negative imaginary a axis. Here, we make 
extensive use of the incomplete Gamma function defined as [see Eq. (6.5.3) of J43|] 



T(a,z) 



dte~ t t a - 1 



(13) 



For a = 0, the quantity T(0, z) as a function of z has a branch cut along the negative real z axis, 
and we assume 



lim lmr(0, — x + ie) 



-7T . 



x>0 



(14) 



which follows from the relationships [43, Eq. (5.1.45)] E\(z) = T(0,z) and [43, Eq. (5.1.7)] 
E 1 (-x±i0) = -Ei(x)=Fivr. 

A unified expansion in terms of special functions - including infinitesimal imaginary parts - 
has also been given in the final version of ||. In this context it is perhaps worth pointing out 
that it is inconsistent with standard notation to use the exponential integral Ei for such a unified 
formula. The exponential integral Ei is usually defined only for real argument. It is defined as a 
Cauchy principal-value integral by Gradshtein and Ryzhik |4JJ [ see Eqs. (8.211.1) and (8.211.2) 
ibid.] as well as by Abramowitz and Stegun [S] [see Eq. (5.1.2) ibid.], and also predominantly in 
the mathematical literature; see, for example, Olver [^6[ [see Eq. (3.07) ibid.]. In contrast to the 
exponential integral Ei, the incomplete Gamma function is defined in the entire complex plane 
with a cut along the negative real axis. 



3 An Important Mathematical Identity 

W. J. Mielniczuk [^7| has outlined a proof of the representation (^) for the real part of the 
effective action. However, his work suffered from a series of unfortunate typographical errors. 
Here, we provide details on the intermediate steps used in our calculation @, and we give, in 
particular, a corrected version of identity (2.8) of [47 1. This corrected version was also used in 
obtaining the results in Eqs. (||) and ( |l0| ) above and in Eqs. (2) — (6) and (18) in J7|. 
The corrected version of identity (2.8) of [47| reads: 
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H > t -^r^ To-^y cotn { - kir ) . (15) 
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In Ramanujan's notebooks |48]| , this identity appears as entry (19.3) on p. 271. A proof is given 
which is based on the repeated application of the well-known identities [see e.g. Eq. (19.2) of @| 
or Eq. (2.4) of @] 

oo 
TTXcot{7rx) = l + 2x 2 ^2l/(x 2 -n 2 ) (16) 

n=l 

and 

oo 
^ycoth(^) = l + 2y 2 ^l/(y 2 +n 2 ), (17) 

n=\ 

respectively, to each one of the factors on the left-hand side of (O), and a skillful reformulation 
of the resulting double sum. 

Here, we will derive an alternative, but fully equivalent formulation of (|l5|): 

1 a 2 - b 2 

ab coth(a z) cot (bz) — —^ — 

z 3 

2abz 2 J y^ coth(A;7ra/6) ^^ coth(A;7rfo/a) [ 

~^~ I _ k I* 2 ~ k 2 TT 2 /b 2 } ~ ^ k [z 2 + k 2 7r 2 /a 2 } I ' (18) 

V. tZ — -L tZ — J. J 

In our derivation of identity ( |l5|) in [|7|], we used the so-called Partial Fraction Theorem 4.4.5 
of J50] . It may be surprising that this identity can be obtained by a straightforward application 
of this basic theorem to the left-hand side of Eq. (|15|), especially in view of the fact that the 
derivation of this result represented a considerable challenge to S. Ramanujan, as it is evident 
from remarks on p. 271 (top of page) of [pq[ . This shows that sometimes interesting new results 
can be obtained by a straightforward application of theorems occurring in standard textbooks. 
We quote here the Partial Fraction Theorem [p0| , Theorem 4.4.5 on pp. 337 - 338)]: 



Theorem 1. Suppose that g is meromorphic with simple poles at a±, ai, 03, • • • with 
< |«i| < I ot2 1 - - ■ and residues b^ at a^, but analytic at 0. Suppose there is an 
increasing sequence R\, R2, R3, ... with linin^oo R n = 00 and simple closed curves 
C n satisfying 

(i): \z\ > R n for all z on C n , 

(ii): There is a constant S with length(C„) < SR n for all n€N, 

(Hi): There is a constant M with \g(z)\ < M for all z on C n and for all n € N. 

Then, 

Here, it should be noted that normally the series on the right-hand side of (^) does not 
converge absolutely, which implies that it must not be rearranged. In fact, the series with the 
terms b m /{z — a m ) may even diverge (compare the remark in the third paragraph on p. 539 
of (5l]]). Thus, the compensatory terms b m /a m are necessary to ensure convergence. 

Equation flil^ ) implies eo ipso - just as the conditions of the partial fraction theorem - that 
the function g{z) which is expanded into "partial fractions" must necessarily be analytic at zero 
argument. In order to analyse the behaviour of the product coth(az) cotipz) with a, b £ R as 
z -» 0, we use [||, Eq. (30:6.2)] 



00 



ooth(*) = i JT w ** ' M<*, (20) 

z 3=0 { J> ' 



and ||, Eq. (34:6.2)] 



■>j 



cot(z) = I-£^^, \ z \<*. (21) 



3 
Here, B m with m € No is a Bernoulli number p3|, Section 23] 



If we now insert the leading terms of (pc|) and (21) into the product coth(az) cot(bz) and use 
B = 1 and B 2 = 1/6 |43|, Eq. (23.1.3)], we obtain: 

coth(az)cot(fe) = — T + ^ - — + O (z 2 ) , z^O. (22) 

This suggests - consistent with the renormalisation of the effective Lagrangian - the definition 
of the function 

1 a 2 - b 2 

f(z) = abcoth(az)cot(bz) ^ , o,kR, (23) 

z A 3 

which corresponds to the left-hand side of Eq. ( |i~8| ) and which satisfies 

f(z) = 0(z 2 ) , z^O. (24) 

For the determination of the poles of f(z), we use the above equations (fig ) and ( |i~7|) which 
can be reformulated as |33, Eq. (30:6.6)] 

00 
coth(z) = >^ — ;t— - , z 7^ ik-K , A;GZ, (25) 



m=—oo 



and m Eq. (34:6.5)] 

CO 

cot(z) = > —z tt—p; , z ^ kir , fceZ. (26) 

v ; ^ z 2 - m 2 7r 2 v y 

m=— oo 

From these expansions and from the definition of f(z) we conclude that f(z) has the simple 
poles 

a k = i — , k£Z\{0}, (27) 

a 

a' k = £jL, fc£Z\{0}. (28) 

Next, we want to determine the corresponding residues 

b k = Res f(z) = lim [(z-a k )f(z)] (29) 

•2— a k z ~ ¥a h 

and 

b' k = Res f(z) = lim [(z-a' k )f(z)] . (30) 



For that purpose, we rewrite ( |25| ) and rt2q) as follows by isolating the terms that contribute 
to the residues at a k and a k , respectively: 

CO , 

t / \ V~^ z / a 
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(z + rmt/b)(z — mir/b) (z + kn/b){z — kir/b) 



V z / 6 + 2z / b (32) 
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With the help of (|3l|), we then obtain: 

h = lim [{z - Ofc) /(z)] 

= lim [(z — ikir/a) aftcoth(az) cot (bz)] 

z^ikn/a 

2bzcot(bz) , , , , . , . 

lim Vt = bcot(ikirb/a). (33) 

z^ikw/a (z + ikir/a) 

If we now use cot(iz) = — icoth(z), we finally obtain 

b k = -ibcoth(k-irb/a), k e Z \ {0} . (34) 

Similarly, we obtain with the help of (|3^): 

ftj. = lim [(z - oi) /(*)] = acot(A;7ra/6) , k G Z \ {0}. (35) 



If ([jl]) is to be used for the derivation of a partial-fraction decomposition for f(z) denned by 
[23|) , it is natural to identify the closed contour C n with the rectangle having the 4 sides X n 



and y\ , where 

4 ±) = K = ^n + %n|x n = (2s - l)X n ,0 < a < l,y„ = ±^n} (36) 

and 

# = {z n = x n + iy n |x n = ±X n , j/ n = (2t - l)Y n , < t < l} . (37) 

Clearly, X n and y n have to be chosen in such a way that the poles a^ = /c7r/6 and a^ = ikir/a 
of /(z) do not lie on this rectangle. 

We are on the safe side if we choose X n and Y n in such a way that they are located in the 
middle between two neighbouring poles of f(z): 

X n = £(n + l/2), Y n = -(n + 1/2), raeN. (38) 

o a 

We now have to show that f(z) is bounded on the rectangle C n according to (iii) of Theorem 1. 
For that purpose, we now use (p2l Eq. (30:11.2)] 



, , , sinh(2x) — i sin(2-u) , . 

cothtx + iy) = -V4 7-^r (39) 

v y> cosh(2x) - cos(2y) v ; 



and ||, Eq. (34:11.2)] 



. sin(2x) -zsinh(2y) 

COt X + IV) = ; : —. -. 4(J) 

v y ' cos(2x) - cosh(2y) v ; 

If s in (36) satisfies s = 1/2, which implies x n = and y n = ±Y n , or t in ( |37| ) satisfies t = 1/2, 
which implies x n = ±X n and y n = 0, we obtain 

, , , . ,, \ , , .,,,x — isin(±2aK,) &sinh(±2frK,) ., . 

coth(±iay n )cot(±t6y n ) = T1 W^Tm tr+o;,v m ' 41 

[1 — cos(±2ar n )J [1 — cosh(±2or n )J 

which remains bounded on C n as n — > oo, or 

coth(±aX n ) cot{±bX n ) = sinh(±2«X w ) sm(±2bX )) 

V ' V 7 [cosh(±2aX n ) - 1] [1 - cos(±26X„)] v ; 

which also remains bounded on C„ as n — ► oo. 



If s in (36) and £ in (j37|) satisfy s,t ^ 1/2, /(z) remains bounded on the rectangles C n 



because of the periodicity and boundedness of cos and sin for real arguments, and we find 

sinh(2ax n ) 
lim coth(a2i n ) = hm — — (43) 

n^oo n^oo COSh(2aX n ) 



and 



lim cot(bz n ) = i lim Sln y,y • (44) 

n->oo n->oo cosh(2oy n ) 

Thus, the constant M in (iii) of Theorem 1 can in the case of f(z) be chosen according to 

M = max(\f(zn)\) , z n G C„ . (45) 



Then, the R n satisfying (i) of Theorem 1 should be chosen according to 

R n = mm(X n ,Y n ), (46) 

and the constant S in (ii) of Theorem 1 should be chosen according to 

S = 4max(6/a,a/6). (47) 



However, the function f(z) defined in (23) does not exactly meet the requirements of Theorem 
1. Instead of a sequence of simple poles oi, 02, 03, • • • with corresponding residues b\, 62, &3j • • • , 
there are actually two sequences o±i, a±2, a±3, ... of simple poles with corresponding residues 
b±i, b±2, b±s, . . . , which implies that we have a partial- fraction decomposition of the following 
kind: 



F W = F(0) + ± {j^ + M + t (t^ + H 

1 4. u m u, m 1 1 xs u '—m u —m J 

m=l v y m=l v ' 



(48) 



If we take into account that the hyperbolic cotangent is odd according to coth(— z) = — coth(z), 
we see from ( p7|) and ( |34"| ) or from ( |28| ) and (35) that the poles and residues of F(z) have to 



satisfy for all k 6 N the (anti-) symmetry relations 

a k = -a_ fe , 6 fc = -6_fc. (49) 

If we also take into account that (p4|) implies 

/(0) = 0, (50) 

we see that ([5J|) is to be replaced by the following partial-fraction decomposition: 



n/ \ v^ I b m b m 2b m I 

F ( z ) = 2^ ) ~z z — + — r 

I 2 O m Z -+- (l m Cl m J 

m=l ' 



(51) 



If we now insert the poles (B7h and ( pq ) and the corresponding residues ([34|) and (|3 
respectively, into (pl|), we obtain the following partial- fraction decomposition: 



1 a 2 - b 2 
ab coth(a z) cot (bz) — —x — — - — 



00 ( 1 12] 

= y^ acoth(/c7ra/6) < ; — — ; — — + - — — > 

-£-J ' ' [z-kir/b z + kir/b k-n/b) 

00 f 1 1 2 1 

- ^ »6coth(fe7r6/a) — —— + —- . (52) 

By putting the expressions in curly brackets on a common denominator, we obtain (|l8|). 

4 Electric— Magnetic Duality 

Electric-magnetic duality, understood as a mutual transformation of electric and magnetic quan- 



tities, has attracted a great deal of attention since Dirac's ideas on magnetic monopoles |54| , |55|| 
were introduced in classical electrodynamics. Whereas electric-magnetic duality can be formu- 
lated as a continuous symmetry of the classical Maxwell equations [EH| (either including both 



types of charges or without any charges), we devote this section to a brief study of discrete 
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duality transformations of quantum effective actions. In particular, we consider two types of 
duality: 

Type I: E -> B, B -> -E, 
Type II: a — > — i6, 6 — ► ia. (53) 

(In a Lorentz frame where -E ||B, the Type II duality implies \B\ — > — i|-E| and |i£| — * i|-B|.) As 
a consequence, the invariants transform as 

Type I: JT — ► —J 7 , Q — » — £/, =>■ 6 — > a, a —>■ —b 
Type II: .F -► F, -» 0. (54) 

Note that the duality of Type II preserves the invariants, so that Maxwell's equations derivable 
from C c \ = —J- are trivially invariant; by contrast, the classical Lagrangian is not invariant 
under Type I, but Maxwell's equations are.g It therefore becomes obvious that only Type I 
affects physical quantities (field strengths), whereas Type II signifies a certain redundancy in 
the parameterisation of physical quantities. 

In fact, this redundancy exists even for a larger class of duality transformations of Type II, 

a— > ±i6 6— >=Fia. (55) 

which preserves the invariants J- and Q. Since the effective Lagrangian of QED is gauge and 
Lorentz invariant, the effective Lagrangian for constant fields is necessarily a function of T and 
Q only, £ e fj = £ e g(JF, Q), implying that the duality transformations of Eq. fl55| ) leave C e Q{T,Q) 
invariant. This statement is not at all tied to perturbation theory, and the duality ( jBlf ) holds 
therefore to all orders in the external fields and to all loop orders. From a diagrammatic perspec- 
tive, each external leg of a diagram contributing to any given loop order corresponds to a field 
strength tensor; and any possible contraction can be expressed in terms of J- and Q, ensuring 
this duality invariance. In some sense, these invariants are therefore more fundamental than a 
and b. 

Let us now study the invariance properties of the effective Lagrangian under consideration. 
From Eq. ([[]), it is obvious that AC is not invariant under Type I, similarly to the Maxwell 
Lagrangian, (even a rotation of the contour in the first quadrant does not help) : in an asymptotic 
expansion, terms of odd order in T, Q flip sign. 

At a first glance, AC indeed seems to be invariant under the transformation ( pq ) as ex- 
pected, since the integrand is invariant.^] Nevertheless, it has recently been argued H[l(| that 
AC is uniquely invariant only under Type II, excluding explicitly the remaining transforma- 
tions of Eq. fl55|). The argument is based on the special function representation ( |T2| ) in terms of 
incomplete Gamma functions. 

In our language, the argument could be rephrased in the following way: Taking the details 
of the integral contour specified by r] and e into account corresponds to a modification a — > 
a exp(i<5), where 5 is a small positive quantity. We discuss the apparent uniqueness of the 
replacement a — > — ib by way of example (see also Appendix |A[). 

When applying the replacement a — > a exp(i 5) to b n in Eq. (|i~2|), the relevant factor in the 
argument of the last Gamma function in b n assumes the form \/a — »• (i/a) exp(— i 6), which has 
complex argument ir/2 — 5. With attention to the fact that the incomplete Gamma function has 
singularities along the negative real axis, it therefore becomes apparent that this argument could 



1 The transformations of a and b under Type I follow from Eqs. (HM. Note that Eq. M is not meant to 
participate in the duality transformations, owing to the nonlinear relation between the two sets of invariants. 

2 The integrand is even invariant under a larger class of duality transformations, including a — > ±ife, b — ► ±i a, 
owing to parity invariance of QED. 
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only be increased by an amount of +ir/2 (not — 7r/2), if we want to stay on the same branch of 
the incomplete Gamma function. This seems to fix uniquely the replacement a — ► — ib because 
in this case, 

i i 2 1 1 

- exp(— i<5) — > — exp(-iJ) — > -- exp(-i£) = —- + 17, (56) 

a b b b 



in agreement with the last term in the definition of c n in Eq. ([i2|). Here, 7 denotes a further 
infinitesimal positive quantity. 

The contradiction between the required general form of the duality (q^) and the seemingly 
unique form of Type II can be resolved by looking at the integral representation: consider the 
representation (|l|) of AC as a function of two complex variables z\ and z 2 : 

oo+ir; 

AC(z 1 ,z 2 ) = — % lim / jL e -(«i 2 -i£)s ZlZ2 co th(e.zis) cot(ez 2 ^ ' ' 



8tt 2 e, V ^0+ J s L 3 ( es ) 2 - ' 

(57) 

The duality of Type II ( |53| ) can be formulated as the identity A£(a,6) = A£(— i6, ia). Because 
a and 6 are real and positive, this identity is valid in a strict sense only if the function AC(zi, z 2 ) 
is one- valued in the relevant range of the complex arguments of z\ and z 2 ■ 

In order to avoid the poles given by the cotangent and hyperbolic cotangent functions, the 
complex arguments of z\ and z 2 in ( |57| ) cannot be varied without restriction if we keep the class 
of allowed contours fixed (recall that the limits of integration in (111) can also be chosen as i 77 and 
oo + ir], respectively). Taking into account the infinitesimal parts e and 77, we are led to conclude 
that the integral representation remains valid for the fixed class of contours in the argument 
range 

-vr/2 < arg(z 1 )<vr/2, 

< arg(z 2 ) < vr (58) 

for z\ and z 2 (observe the fine difference between the < and < signs!). The restriction given in ( |58| ) 
identifies (as a function of Z\ and z 2 ) a "physical" or "causal" branch of the effective action for a 
given contour. Among the four different replacements AC(a, b) — > A£(±i6, ±ia) or AC(a, b) — > 
A£(±i6, =Fia), h is only the replacement AC(a, b) — » A£(— i&, ia) which respects the restriction 



set by (58) for a fixed contour. Therefore, the duality of Type II is only "unique" in connection 
with a precisely specified integral representation which does not allow for an unrestricted analytic 
continuation of its arguments. In other words, the seeming "uniqueness" is simply a shortcoming 
of the integral representation (and also of the special-function representation being identical to 
the former). 

In order to achieve full invariance under the dualities ( |55[ ) , we have to allow for the fact that 
the contour also has to be readjusted (or certain poles are allowed to be crossed without picking 
up their contribution). From a different perspective, it is only natural to perform the duality 
transformation first, and then specify the details of the contour. This is perfectly justified, since 
the particular choice of the contour is not a result of the calculation (of the fermion determinant 
in this case), but rather an additional piece of information that has to be inserted afterwards 
in order to define the integral. This information arises, of course, from physical requirements: 
the ±ie prescription is dictated by causality, and the shift by rj can, e.g., be fixed by requiring 
that the pair-production probability related to the imaginary part of AC is a number between 
and 1. The resulting integral representation will always lead to the same special- function 
representation (|l^). 
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From another point of view, the special-function representation and the integral representa- 
tion with fixed contour remove a part of the above-mentioned redundancy in the parameterisa- 
tion of the field strength invariants, which is generally present in the effective Lagrangian. 

At this point, let us remark that dualities of Type II or (^) can be very useful from a 
technical viewpoint, although they have no physical meaning: for instance, from the effective 
Lagrangian for a purely magnetic field, one can extract information about the electric case with 
the aid of the substitution, \B\ — > — i\E\ [63]. Moreover, standard-model calculations in constant 



electromagnetic fields can always be checked by testing their dual invariance; e.g. in [57, 5^], this 
dual invariance is visible in neutrino amplitudes in electromagnetic fields. Of course, to be on 
the safe side, the duality of Type II suffices for such a check in order to avoid problems of the 
kind mentioned before. 

The question of electric-magnetic duality becomes even more interesting for systems which 
are characterised by additional Lorentz covariant quantities. For example, let us consider QED 
with constant fields in a heat bath; the latter involves an additional Lorentz vector Uu, the heat- 
bath four- velocity, allowing for one further gauge and Lorentz invariant quantity [cf. Eq. (3.143) 
of@]: 

£ = ( Ufl F^) 2 = u^ v uPF pv . (59) 

Under duality of Type II, £ is trivially invariant, since it is linearly independent of J- and Q, 
and thereby independent of a and b. In fact, the known QED effective actions, depending on a, 
b and £ at finite temperature, show this invariance under Type II Q. 

Under duality of Type I, £ transforms according to £ — ► £ + 2J 7 ; and similarly to the zero- 
temperature case, the finite-temperature effective action is generally not invariant under Type I. 
Nevertheless, the dominant low-temperature contribution arising at two-loop order astonishingly 
exhibits an invariance under Type I. This might be related to the fact that only transversal 
thermal fluctuations of the photon give rise to this contribution. 

5 Acceleration of Convergence 



From the asymptotic expansion of the incomplete Gamma function as z — ► oo 43, Eq. (6.5.32)], 
e z T(a,z) = z a - 1 



l + ^ + (a -^- 2) +0 (,- 



z^ 



|argz| <3vr/2, (60) 



we obtain 

e nz T(0,nz) 



n n 2 z n 3 z 2 



+ 0(n~ 4 ) , n->oo. (61) 



Since the cotangents in (|l2| ) rapidly approach one as n —>■ oo, we can conclude that terms in ([12] ) 
are essentially of order 0(n~ 2 ) as n — > oo. A more detailed analysis of the large-order behaviour 



of the terms in (12) can be found in Appendices A and B of J/J. 

An 0(n -2 ) behaviour of the terms also occurs in the Dirichlet series £(s) = Y^n=\ n ~ s f° r the 
Riemann zeta function with s = 2. This is a very discomforting observation since the Dirichlet 
series for ("(2) converges quite slowly - it can be shown that increasing the number of terms of 
its partial sum by a factor of 10 improves the accuracy by a single decimal digit only - and in 
the literature on convergence acceleration it is one of the standard test systems for checking the 
capability of a transformation in the case of a slowly convergent monotonic series. 

Thus, the series expansions @ and dl2| ) represent typical cases of monotonic series which 
converge so slowly that a straightforward evaluation by adding up one term after another is 
computationally very unattractive. Of course, the acceleration of the convergence of series of 
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that kind has been studied quite extensively in the literature and many techniques are known 
to improve the efficiency of numerical computations (see for example [59 plfl and references 
therein). 

Nevertheless, one should not forget that the acceleration of the convergence of a slowly 
convergent monotonic series can be a very challenging problem. Moreover, numerical instabilities 
due to rounding errors are likely to occur. A sequence transformation accelerates convergence 
by extracting and utilising information from a finite set of input data on the index-dependence 
of the truncation errors. This is normally accomplished by forming higher weighted differences. 
If the input data are the partial sums of a strictly alternating series, the formation of higher 
weighted differences is a remarkably stable process, but if the input data all have the same sign, 
numerical instabilities are quite likely. Thus, if the sequence to be transformed are the partial 
sums of a slowly convergent monotonic series, numerical instabilities are to be expected, and 
most convergence acceleration methods are not able to obtain transformation results that are 
close to machine accuracy. 

In some cases, these instability problems can be overcome with the help of a condensation 
transformation attributable to Van Wijngaarden, which converts input data having the same 
sign to the partial sums of an alternating series, whose convergence can be accelerated more 



effectively (compare, for instance, Appendix A of [42]). The condensation transformation was 
first mentioned in |6^, pp. 126 - 127] and only later published by Van Wijngaarden Q. It was 



used by Daniel [64] in combination with the Euler transformation, and recently, it was rederived 
by Pelzl and King |35| . Since the transformation of a strictly alternating series by means of non- 
linear sequence transformations is a stable process, it was in this way possible to evaluate special 
functions that are defined by extremely slowly convergent monotonic series, not only relatively 
efficiently but also close to machine accuracy [42], or to perform extensive quantum electrody- 
namical calculations |4y,^l[. Unfortunately, the use of this "combined nonlinear-condensation 
transformation" (CNCT) p2| is not always possible: The conversion of a monotonic to an alter- 
nating series requires that terms of the input series with large indices can be computed. 

However, this CNCT is well suited for the acceleration of the convergence of the series 
expansions @ and fll2|) considered in this article. 

The method consists in first rewriting the slowly convergent monotonic input series 
Sfclo r (^) ^° an alternating series. In the second step, the convergence of the alternating 
series is accelerated via a suitable nonlinear sequence transformation. In our case, we have 
r{k) = [afc+i + ^fc+i]> with a^+i and c4+i given by Eq. (0). The partial sums 



ov, 



E r ( fc ) 



(62) 



fc=0 



increase monotonically, i.e. <7 n +i > o n for all n = 0, . . . , oo if all terms satisfy r(n) > 0. Let us 
further assume that the sequence of the partial sums {a n }^ =0 converges to some limit a = o"oo 
as n — > oo. Following Van Wijngaarden [63], we transform the original series into an alternating 
series Y1T=q(~ 1) j A? according to 



£>(*) 



k=0 



A.i) 



£(-i) j a,, 

oo 
fc=0 

2 fe r(2 fc (i + l)-l). 



(63) 

(64) 
(65) 



The terms Aj defined in Eq. (p4|) are all positive if the terms r(k) of the original series are 
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Table 1: Evaluation of the real (dispersive) part of the QED effective Lagrangian 
for \E\ = 30 -E^ and |B| = 30i? cr by evaluating the special function representation 
(g) using a combination |4^] of the Van Wijngaarden condensation transformation 



defined in Eqs. ( p3|) - ( |65| ) and the nonlinear delta transformation [59]. The S n are 
defined in (|6q), and the nonlinear delta transform S n (l, So) is defined in Eq. (8.4-4) 
of |J59]| . The result is given in units of m 4 . The apparent convergence of the delta 
transforms in the second column is indicated by underlining. After 21 iterations, the 
transforms <5.i 0) (l,S ) have stabilised to the 15-figure result 8.393 398 582 100 617. 



n 


^>n 


<4 U) (1,S ) 





11.834 710 587 368 


11.834 710 587 368 


1 


6.388 353 476 842 


8.463 873 830 587 


2 


9.741 830 922 440 


8.384 963 963 657 


3 


7.413 293 009 436 


8.393 553 703 382 


4 


9.141 317 648 944 


8.393 398 289 155 


5 


7.803 272 984 326 


8.393 399 592 701 


6 


8.870 519 831 631 


8.393 398 337 299 


7 


8.000 385 721 330 


8.393 398 666 561 


8 


8.721 936 912 365 


8.393 398 594 227 


9 


8.115 447 030 287 


8.393 398 583 148 


10 


8.630 906 154 429 


8.393 398 580 473 


11 


8.188 731 860 505 


8.393 398 582 040 


12 


8.571 048 192 682 


8.393 398 582 143 


13 


8.238 223 421 331 


8.393 398 582 104 


14 


8.529 698 152 551 


8.393 398 582 097 


15 


8.273 084 942 867 


8.393 398 582 099 


16 


8.500 073 809 529 


8.393 398 582 099 



all positive. The quantities Aj are commonly referred to as the condensed sums, and the series 
^°^ (— l) J Aj is referred to as the Van Wijngaarden transformed series. 

The transformation from a monotonic series to a strictly alternating series according to 



Eqs. (p3[), (64) and ( |6q ) is essentially a reordering of the terms r(k) of the original series. We 



define the partial sums 



S n = £(-1)^ (66) 

3=0 



of the Van Wijngaarden transformed series. As illustrated in Table 1 of [[43], the S n with n > 
reproduces the partial sum a n , Eq. (|62|), which contains the first n + 1 terms of the original 
series. Formal proofs of the correctness of this rearrangement can be found in Ref. |54j or in the 
Appendix of Ref. |65| . 



The series (64) for the terms of the Van Wijngaarden transformed series can be rewritten as 
follows: 

Aj = r(j) + 2r(2j + 1) + 4r(4j + 3) + . . . . (67) 

Since the terms r(k) of the original series have by assumption the same sign, we immediately 
observe 

|A.j| > |t-C/)I- 
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Consequently, the Van Wijngaarden transformation, given by Eqs. (pq), ( |64| ) and (p5|), does not 
lead to an alternating series whose terms decay more rapidly in magnitude than the terms of the 
original monotonic series. However, an acceleration of convergence may be achieved if the partial 
sums ( |66| ) of the Van Wijngaarden transformed series are used as input data in a convergence 
acceleration process, and - as, for example, discussed in Appendix A of |42|]) - the convergence 
of alternating series can be accelerated much more effectively than the convergence of monotonic 
series. For the acceleration of convergence we use the delta transformation given in Eq. (8.4-4) 
of pjj, which was found to be a very effective accelerator for Van Wijngaarden transformed 



series |42|| . 



This will now be demonstrated explicitly. From a consideration of the expressions (||) and 
( |i~2| ) it is obvious that the computationally most demanding special cases will be encountered for 
large fields; in these cases, many terms of the representation (|9]) are needed in order to achieve 
convergence [accordingly, for strong fields we encounter problematic oscillations in the integrand 
of Eq. (U)]. We consider only one example here - Table || -, which describes the evaluation 
of the dispersive (real) part of the effective Lagrangian at field strength \E\ = b = 30E cr 
and \B\ = a = 30B cr . This does not preclude the possibility that other efficient calculational 
methods exist for the evaluation of (|9]). However, we stress here that a suitable acceleration 
of the convergence of the special function representations (g) and (^) removes the principal 
numerical difficulty associated with the slow convergence of the series at large field strength. In 
our example - see Table [l] -, the highest index encountered in the calculation is h = 37748736, 
the total number of evaluations of terms [a n + d n ] is 405. The ratio of these two numbers is 
roughly 93000, which corresponds to an acceleration of the calculation by roughly five orders of 
magnitude. 



The CNCT transforms the slowly convergent series expansions @ and (12) into the rapidly 
converging sequence of the delta transforms <% (l, So) (see Table |l|). In general, no direct in- 
terpretation is available for the delta transforms []59|, just as much as Pade approximants p7[ ] 
lack a direct physical interpretation. At best, the delta transforms can be viewed as the analytic 
continuations ("extrapolations") of the partial sums (p6|) of the Van Wijngaarden transformed 
alternating series S n to n — > oc. 

6 Conclusion 

We have investigated questions related to the representation of the quantum electro dynamic 
(QED) effective Lagrangian and its numerical evaluation. In Sec. 0, we recall our previous 
results given in Qb] for special function representations of the effective Lagrangian, and we 
clarify the mathematical notation used in the special function representations (B), (llO) and 
(|12|) . The representation (12) unifies the real and imaginary parts. 



In Sec. ||, we present a detailed description of the proof of a certain mathematical identity 
( J15f) used in our investigations @,g. In Sec. ||, we discuss the question as to whether the QED 
effective Lagrangian is invariant under certain types of electric-magnetic duality. We conclude 
that gauge and Lorentz symmetry guarantees invariance under a general class of dualities (|55| ) 
to all loop orders; but discrete representations of the effective Lagrangian may not realize this 
general dual invariance in a strict sense, so that only a smaller duality subgroup (Type II) 
remains. 

In Sec. pi, we show that the convergence over the summation index n of the special function 
representation (||) can be accelerated efficiently by the application of the CNCT transforma- 
tion p2| |. In this way, the computing time is reduced by several orders of magnitude. Based 
on the results of the current paper, we expect to carry out detailed studies related to various 
projected and ongoing experiments and astrophysical phenomena |l, 12-14,16] involving strong 
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static-field conditions (or fields with frequencies that are small as compared to the electron 
Compton wavelength). 
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A Model Example 

Consider as a model example the integral 

[■oo+ir, / 1 1 1 1 \ 

M(a, b) = di exp(-t) : — + : — + — + — , (69) 

K,J J ir) FV ' \l-iat 1 + iat l-bt l + btj' v ; 

where a and b are a priori real variables, but can be also be generalised to the complex case 
a — ► zi, b — ► Z2- The integrand is invariant under the four different "duality transformations" 
a — ► ±i&, b — » ±ia (=Fia), but the integral is not. Owing to the poles of the integrand, the 
imaginary parts change sign when the complex arguments a and b are varied such that one of 
the poles of the integrand is crossed. When considered as a function of complex arguments z\ 
and Z2, the function M(z\, Z2) has branch cuts along the positive and negative imaginary z\ axis 
and along the positive and negative real Z2 axis. Furthermore, by inspection of the integrand we 
conclude that the argument ranges for z\ and Z2 are 

-tt/2 < arg(3i)<7r/2, 

< arg(z 2 ) < vr . (70) 

From these relations it becomes clear that if we want to stay on the principal branch in the 
complex Z\,Z2 plane, then we have to modify the arguments of a and b in accord with the 



restrictions given by (|7C]). This singles out the duality of Type II 

a — > — ib , b — > ia . (71) 

We conclude the discussion of the model example by pointing out that it can be expressed as 

exp (D r (°^)- cxp (4) r ( '4 



M(a, b) = lim 

e^0+ I a 



<X!> « I) r ( '^) - exp (-0 r (°'-l +ie 



(72) 



1 
+ b 

There is a certain analogy to the special function representation fll2|). 
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